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Abstract. We have studied the two-body production and decay of a heavy, right- 
handed neutrino to two hght states using the full Boltzmann equation instead of the 
usual integrated Boltzmann equation which assumes kinetic equilibrium of all species. 
Decays and inverse decays are inefhcient for thermalising the distribution function of 
the heavy neutrino and in some parameter ranges there can be very large deviations 
from kinetic equilibrium. This leads to substantial numerical differences between the 
two approaches. Furthermore we study the impact of this difference on the lepton 
asymmetry production during leptogenesis and find that in the strong washout regime 
the final asymmetry is changed by 15-30% when the full Boltzmann equation is used. 
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1. Introduction 

Leptogenesis is perhaps the most attractive model for generating the matter-antimatter 
asymmetry in our Universe [1] after inflation. The process generates a lepton asymmetry 
via the production and subsequent decay of a heavy Majorana neutrino. This lepton 
asymmetry is partially converted to a baryon asymmetry via sphaleron processes [2] 
which break both B and L, but conserve B — L. 

The fraction oi B — L that ends up in B by the sphaleron processes is given by 
O'sph — 28/79 giving ub = cisphnB-L — —dsphi^L where the ub-l is created by leptogenesis 
which gives a non zero value of L (see for instance [3]). 

In its simplest form (within the context of the see-saw model) leptogenesis consists 
of adding three heavy right-handed neutrinos to the standard model. In the hierarchical 
limit one of these right handed neutrinos, vr^, is much lighter than the other two and 
the leptogenesis mechanism consists essentially of the process 



where is a light, left-handed neutrino and (f) the Higgs. All light flavours behave 
identically, so you can think of the other light flavours as included as a factor of 3 in 
the decay rate. Because of loop corrections there can be CP violation in the decay, 
usually quantified by the asymmetry parameter e. In the following we denote the heavy 
neutrino by R, the hght by L, and the Higgs by H. 

This model has been studied extensively in the literature, including deviations 
from the hierarchical limit, thermal corrections etc [3-17]. However, all studies have 
used the integrated Boltzmann equation to follow the evolution of the heavy neutrino 
number density and the lepton asymmetry. This approach assumes Maxwcll-Boltzmann 
statistics for all particles as well as kinetic equilibrium for the heavy species. This 
assumption is normally justified in freeze-out calculations where elastic scattering is 
assumed to be much faster than inelastic reactions. However, in the present context, 
kinetic equilibrium in the heavy species would have to be maintained by the decays and 
inverse decays alone. Therefore it is not obvious that the integrated Boltzmann equation 
is always a good approximation. Furthermore 1 <-> 2 processes are generally inefficient 
for thermalization compared with 2^2 processes. For 2^2 deviations from kinetic 
equilibrium are always of order 20% or less [18], but for 1 ^ 2 processes they can be 
very large (see for instance Refs. [19-22] for a case where deviation from equilibrium is 
extremely large). 

In this paper we investigate how the use of the full Boltzmann equation affects the 
final lepton asymmetry in a simplified model with only decays and inverse decays and 
resonant scattering. We will return to the point of resonant scattering in due course. We 
find that when T ~ tur the difference can be very large. However, at small temperature 
where the inverse decay dominates the difference decreases in magnitude to about 20%. 




(1) 



Decay of heavy Majorana neutrinos using the full Boltzmann equation including its implications for leptogent 



2. The Boltzmann equation 

Here we study only the two-body decay of a heavy right-handed neutrino to a light 
neutrino plus a Higgs. We do not include thermal corrections to the particle masses, so 
that for instance the process H — > LR is not kinematically possible. We assume that 
the asymmetry, represented by 
r - f 

is small, so that when we calculate anything with R, we can assume identical 
distributions of L and L. We consider only initial zero abundance of R. We use 
only single particle distribution functions, in which case the Boltzmann equation for 
the heavy species can be written as [4, 23] 
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and for the light neutrino it is 
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with a similar equation for L. The interesting Boltzmann equations for the present 
purpose are those for R and for B — L. H and L, L have gauge interactions which are 
very fast. This means that H can be described by a distribution in chemical equilibrium, 
and that L, L can be described as distributions in kinetic equilibrium, 

/^=(l + e^«/^)-^ (5) 

f^^{l + e^P^-^y^)-' (6) 

f^ ^ (I + e(PL+^y^)-\ (7) 

with ii/T = 3{nL - ni)/T^ + 0{{ix/Tf). 

Using CPT-invariance, following the idea of [23] , we find 



\Mr_,hl\ 
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This simplifies the Boltzmann equations. 

For the heavy neutrino the Boltzmann equation can be simplified to [20] 



dfR pjj^fn 
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dpH [fHfL{l - fR) - fR{l - /l)(1 - fn)] , 
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where Ftot = ^r-^lh + ^r-^lh = T + F is the total rest frame decay rate. The 
corresponding equation for L — L \ is 

— di — dp 

"^""^ C dEn[(fH + fR)F- 



2ElPl J^+pl 



+ e( - 2(1 + fH)fR + Ur - Ml - 2fR))F^ 



(11) 



to first order in e and with 



F^^fL + fL = + ^((Z^/^)') (12) 



2ef/^ 



(14) 

Since wc have so far only included 2^1 processes, Eq. 11 suffers from the well- 
known problem of lepton asymmetry generation even in equilibrium [23]. To remedy 
this problem the resonant part of the LE. <-> LE. must be included. To lowest order in 
e this amounts to adding the term 2e(l — /R)/ff-F"^ in Eq. 11 [8], so that the final form 
of the equation for the lepton asymmetry is 

^Ul - fl) ^^ difL-fi) .... 
dt d^ 

- fl dEn UfH + fR)(F- + eF+) - 26/^(1 + fn) 

This equation docs not exhibit any lepton asymmetry generating behaviour in 
thermal equilibrium because [{fn + /i?)(-F" + eF^) — 26/^^(1 + fn)] = explicitly. 

To first order in /x/T equations Eq. 10 and 16 can be easily integrated numerically 
for given values of tur and F. If Maxwell-Boltzmann statistics and kinetic equilibrium 
for it! are assumed the equations can be further simplified. The integrated Boltzmann 
equations are then 



nR + 3HnR = - (F) (n^ - n^^eq) (16) 
riL-L + ^HnL_i = e (F) {ur - nR,eq) + ^/Ci (m^/r) ^ , (17) 

where z = ^ and (F) = {T + T)ICi{z) / IC2{z) is the thermally averaged total decay rate. 

These equations are the ones normally used in leptogenesis calculations for the 
simplest case of one massive and one light neutrino. However, compared with 
Eqs. (10-16) they involve the approximation of assuming Maxwell-Boltzmann statistics 
and kinetic equilibrium. Particularly for the case where T/H{T — tur) ~ 1 this 
approximation is not necessarily good. 

I Note that wc ignore the sphaleron L to B conversion during the decay process because we track only 
L — L. Including it would have only a modest effect on our numerical results. 
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When decays and inverse decays are included (as well as gauge interactions of L and 
H which establish kinetic equilibrium for L, L, and H individually, but which conserve 
L — L) there are essentially only two interesting parameters. All terms that contribute 
to the development of densities are proportional to either H or F. This competition can 
be parameterised by the single parameter K = T/H(T = niR), the decay parameter. 
The other important parameter is the net decay asymmetry e. 

In the present study we assume that the heavy neutrinos start with zero abundance, 
i.e. that they are not equilibrated at high temperatures through other interactions. We 
refer the reader to [8] for a thorough discussion of the various possibihties for initial R 
abundance. 

In Figs. 1-3 we show a calculation of nji and the net asymmetry for different values 
of K. At high temperatures the equilibration rate of R is much higher when the full 
Boltzmann equation is used. The main reason for this can be seen in Fig. 4. When 
z < 1 the low momentum states of R are populated very efficiently, while there is 
little population of high momentum modes. This in turn means that the inverse decay 
HL — > i? is suppressed relative to the decay process R HL. 

To quantify this, we want to define a momentum dependant inverse decay rate. 

The efficiency with which R states with momentum pn are produced from the 
background (ignoring decays and Pauli Blocking) is then given by 
T{pr) 1 dfn mnV r{Eu+Pn)/2 



/ dpHfnfL. (19) 
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This function is sharply peaked at low momentum with the following limiting behaviour 

F 
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(20) 



This is a steeply increasing function of z and at low temperatures when T « ruR 
the distribution of R is in complete equilibrium and << 1 so that the Boltzmann 
approximation is valid. Therefore the curves for nR/ng^ approach each other in Figs. 1- 
3. Note, however, that the highest pr modes are always out of equilibrium because of 
the asymptotic 1/pr term. In fact for given values of z and K there will be a limiting 
value of Pr above which the distribution will be out of equilibrium. In Fig. 5 we show 
the value oi pr/T for which V{pr)/H = 1 as a function of z, for the specific case of 
K — 1. At high temperatures no momentum states are in equilibrium, whereas at 
low temperatures only progressively higher momentum states are out of equilibrium. 
This figure corresponds well to what is shown in Fig. 4. For z = 0.2 no states are in 
equilibrium, for 2; = 1 states above Pr/T ~ 1 are not yet in equilibrium, and for z = 5 
almost all states have been populated to their equilibrium value. 
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Figure 1. Upper panel: The evolution of rijiln^^^ for the full (solid line) and integrated 
(dotted line) Boltzmann equations. Bottom panel: The evolution of I'^l-lI/^'^ fo'" the 
same cases. The calculation is for X = 0.1 and e = 10^^. 



As can be seen in Figs. 1-3 the net asymmetry also grows much more rapidly 
initially when the full Boltzmann equation is used, and the change of sign occurs at 
higher temperature. However, at low temperatures the differences are smaller. 

In Fig. 6 we show the ratio of the final asymmetry in the two cases 





- full 







As can be seen from Fig. 6 there can be a significant difference in the asymmetry 
between the two approaches. For very high values of K the difference is small because 
the distribution is kept very close to kinetic equilibrium (for a Fermi-Dirac distribution) 
throughout the decay. Note, however, that the difference never goes to zero. The 
reason is that the thermal decay width is changed relative to the Maxwell-Boltzmann 
approximation when Pauli blocking and stimulated emission factors are included (as 
also noted by [8]). For smaller values of K the difference increases and can be as large 
as 50%. However, it should be noted that our framework will break down for small 
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Figure 2. Upper panel: The evolution of rifl/ncq for the full (solid line) and integrated 
(dotted line) Boltzmann equations. Bottom panel: The evolution of I'^l-lI/^'^ fo'" the 
same cases. The calculation is iai K = \ and e = 10^^. 



values of K because we have not included scattering. 

Another important point to note is that r^fuu/^^mt does not depend on e because the 
deviation from kinetic equilibrium is governed completely by the total decay rate (i.e. 
by K)^ not by the asymmetry, as long as e << 1. 

We are only interested in the final value of the asymmetry since, in order to 
maximise the baryon asymmetry for any given value of e, we want leptogenesis to end 
before the conversion through sphalerons end. 



3. Discussion 



We have solved the full Boltzmann equation for decays and inverse decays to follow the 
generation of lepton asymmetry during leptogenesis. When decays are semi-relativistic 
the difference between using the full Boltzmann equation and the standard integrated 
Boltzmann equation can be very large. However, at low temperature where washout 
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Figure 3. Upper panel: The evolution of nn/ricq for the full (solid line) and integrated 
(dotted line) Boltzmann equations. Bottom panel: The evolution of In^^i/T^l for the 
same cases. The calculation is for K = 100 and e — 10~^. 



dominates the difference is relatively modest. 

The difference at K >1 is less than about 30% between the two approaches. For 
smaller K the difference can be larger. However, this is the regime where 1^2 processes 
do not dominate, i.e. outside the regime in which our framework is valid. 

The conclusion is that the Boltzmann approximation yields results in the strong 
washout regime which are accurate to at least 30%. For K > 5 the difference is not 
larger than 15%. 
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is for K = 1. 
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